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ABSTRACT 


The  plane  strain  problem  of  a  step  load  moving  with  uniform  superseismic 
velocity  V  >  Cp  on  the  surface  of  a  half -space  is  considered  for  an  elastic- 
plastic  material  obeying  the  von  Mises  yield  condition. 

Using  dimensional  analysis  the  governing  quasi-linear  partial  differential 
equations  are  converted  into  ordinary  nonlinear  differential  ones  which  are  solved 
numerically  using  a  digital  computer.  To  overcome  computing  difficulties  asymp¬ 
totic  solutions  are  derived  in  the  vicinity  of  a  singular  point  of  the  differ¬ 
ential  equations. 

Numerical  results  are  presented  for  a  range  of  selected  values  of  the 

po 

significant  ncndimensional  parameters,  i.«r.  of  the  surface  load  —  ,  of  Poisson's 

V 

ratio  V  and  of  the  velocity  ratio  —  . 
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LIST  OF  SYMBOLS*  ^ 

Functions  defined  by  Eqs.  (A-8)  to  (A-15). 

Functions  defined  by  Eqs.  (III-l4)-(III-l6)  and  (111-26). 

Velocity  of  propagation  of  elastic  P-waves,  S-vaves  and  of 
inelastic  shock  fronts,  respectively. 

Plastic  potential,  Eq.  (il-l). 

Shear  modulus. 

Invariants  of  stress. 

Yield  stress  in  shear,  Eq.  (II-l). 

Bulk  modulus . 

Function  related  to  inelastic  behavior,  Eq.  (11-31). 

Surface  pressure. 

Intensities  of  step  load  surface  pressure. 

Principal  stress  deviators. 

Stress  deviators  with  respect  to  axes  x,  y,  etc. 

Time. 

Components  of  particle  velocity  in  x  and  y  directions, 
respectively. 

Velocity. 

Cartesian  coordinates,  Fig.  1. 

Jfondimensional  expression. 

Values  of  X  at  P-  and  S-fronts,  respectively. 

Hondimensional  stress  variable. 

Angle  between  the  directions  of  s1  and  of  the  position  ray  of 
an  element.  Fig.  3* 


Other  symbols,  which  are  used  in  one  location  only,  are  defined  as  they  occur 
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A  ■  3  -  3 
Ac,  Av,  At,  etc. 

«  *  q>  -  9 

*ij  ;  *i<5 
T)  ■  Y  -  £ 

6 

X  >  0 
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5  ■ x  z.  yt. 
y 

p 

°u  *  4u 

*1  '  a2  '  °3 
T 

<? 

«pp  >  9s  >9 

*1  '  *2  '  *3  '  ^ 


Small  quantity  for  purposes  of  asymptotic  expansion. 
Increments  of  ff,  v>  T,  etc.  at  a  front. 

Small  quantity  for  purposes  of  asymptotic  expansion. 

Strains,  strain-rates. 

Small  quantity  for  pus-poses  of  asymptotic  expansion. 

Angle  defining  direction  of  the  principal  stress  ,  Fig.  3» 
Function  related  to  inelastic  behavior,  Eq.  (II-8). 

Poisson's  ratio. 

Variable  defined  by  Eq.  (II-l6). 

Mass  density  of  medium. 

Stresses,  stress  rates. 

Principal  stresses. 

Shear  stress. 

Position  angle  of  element,  Fig.  4. 

Position  of  the  elastic  P-  and  S-  and  of  the  inelastic  shock 
fronts,  respectively. 

limits  of  inelastic  regions. 
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I  INTRODUCTION. 

The  two  dimensional  problem  of  the  effect  of  a  pressure  pulse  p(x  -  Vt) 
progressing  with  the  velocity  V  on  the  surface  of  an  elastic  half -space,  Fig.  1, 
has  been  treated  by  Cole  and  Huth  [1J  for  a  line  load  and,  by  superposition,  may 
be  found  for  any  other  distribution  p(x  -  Vt).  Miles  [2J  has  considered  the  three 
dimensional  problem  of  loads  with  axially  symmetric  distribution  p(r,t)  over  an 
expanding  circular  area  on  the  surface,  Fig.  2.  He  has  demonstrated  that  the  plane 
problem  [1J  contains  the  asymptotic  solution  for  the  three  dimensional  problem  [2J 
in  the  region  near  the  wave  front.  The  actual  solution  of  the  three  dimensional 
problem  would  require  a  great  numerical  effort  which  can  be  avoided,  by  using  the 
solution  of  the  p3.ane  problem  to  estimate  the  effect  of  circularly  expanding 
surface  loads. 

Real  materials  can  not  be  expected  to  be  elastic,  and  solutions  of  the  three 
dimensional  problem,  Fig.  2,  for  dissipative  materials  are  hopelessly  complex. 
However,  estimates  for  the  three  dimensional  case  can  be  made  from  generalizations 
of  the  problem  treated  in  [1J  for  dissipative  materials.  This  has  been  done  for 
linearly  viscoelastic  materials  by  Sackman  [3'J,  and  Workman  and  Bleich  [4j,  in  the 
superseismic  and  subseismic  ranges,  respectively.  Based  on  a  formal  solution  in 
[5]  the  effect  of  a  step  load  moving  with  superseismic  velocity  is  determined  in 
the  present  report  for  an  elastic-plastic  material  obeying  the  von  Mises  yield 
condition.  The  identical  problem,  for  a  yield  condition  suitable  for  materials 
with  internal  slip  subject  to  Coulomb  friction,  is  concurrently  being  treated  for 
publication  elsewhere  [6 J . 


The  yield  mechanism  in  the  medium  makes  the  problem  nonlinear,  such  that 
superposition  is  not  permitted  and  each  pressure  distribution  p(x  •>  Vt)  poses  a 
separate  problem.  The  present  paper  treats  the  case  of  a  progressing  step  load 


RHN 
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p(x  -  Vt)  »  p  H\Vt  -  x).  An  approach  permitting  an  approximate  oolution  of  the 
important,  but  very  complex  case  of  a  decaying  surface  pressure  is  discussed  in 
the  Conclusions. 

The  formulation  of  the  problem  in  [5]  furnishes  a  set  of  simultaneous, 
ordinary  nonlinear  differential  equations  which  are  solved  numerically  by  a 
Runge-Kutta  forward  integration  scheme  [73  utilizing  an  IBM  7090*  The  integration 
encounters  numerical  difficulties  near  singularities  of  the  system  of  differential 
equation.  To  treat  such  situations  asymptotic  expansions  were  employed,  leading 
to  approximate  differential  equations  which  could  then  be  integrated  numerically. 

Numerical  results  for  all  combinations  of  the  nondimensional  parameters 

v  =  0,  t  ,  r-  ,  0.35  and  ~  «  1.25,  1*5>  2.0,  if.0,  are  given  in  Tables  l-l6  in  each 
o  4  Cp  p 

case  for  five  values  of  the  surface  load  .  The  coverage  of  results  is  extensive 
enough  to  permit  interpolation. 
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11  FORMULATION  OF  THE  BASIC  AQUATIONS. 

Figure  3  indicates  the  half -space  and  a  system  of  Cartesian  coordinates.  The 

x-axie  is  in  the  direction  of  motion  of  the  step  load,  the  y-  and  z-axes  are  normal 

to  the  surface  in  and  out  of  the  plane  of  the  figure,  respectively.  The  analysis 

considers  the  case  of  plane  strain,  e  *  0,  when  the  velocity  V  of  the  step  load 

z 

is  superseismic.  It  is  known  [8]  that  the  largest  wave  velocity  in  an  elastic- 
plastic  material  is  Cp  ,  the  velocity  of  elastic  P-vaves  in  the  material,  so 
that  the  terra  superseismic  means  V  >  cp  .  Throughout  the  analysis  it  is  assumed 
that  the  strains  and  velocities  are  small,  so  that  their  higher  powers  may  be 
neglected  in  comparison  to  linear  terms. 


To  describe  the  behavior  of  the  elastic-plastic  material  the  plastic  potential 
is  introduced 


»  p 

f  ■  j2  - 

s 

where  Jp  is  the  invariant 

J2  -  I  Sij 

and  the  value  k  >  0  is  the  yield  stress  in  shear. 


(II-l) 


(II-2) 


The  behavior  of  an  element  of  the  material  is  then  defined  by  the  statements 
which  follow. 

1.  The  value  of  the  function  F  may  never  be  positive 

7  <  0  (II-3) 

2.  If,  in  an  element  of  the  material  at  a  given  instant, 

7  <  0  (II-4) 

the  rates  of  change  in  stress  and  strain  are  related  by  the  conventional 


elastic  relations 


However,  if  the  yield  condition 


F  *  0  (II-5) 

is  satisfied  three  possibilities  exist:  a)  in  the  next  instant  of 
time  the  material  may  be  in  a  state  of  plastic  deformation;  b)  it 
may  be  in  a  state  of  elastic  unloading;  c )  it  may  be  in  a  neutral 
state . 


a)  If  the  material  is  in  a  state  of  plastic  deformation 


f  =  0 


(II-6) 


the  total  strain  rate  will  be  the  sum  of  an  elastic  and  a 
plastic  portion 


*  •£) 

•id  ’  •id 


+  « 


id 


(II-7) 


•E 

where  is  obtained  from  the  conventional  elastic  relations, 
.while 


X,  which  must  be  positive 


X  >  0 


(II-8) 


(II -9) 


is  an  a  priori  unknown  function  of  space  and  time.  It  is  to 
be  found  as  part  of  the  solution  of  the  problem. 

b )  In  case  of  elastic  unloading  F  <  0  holds  and  the  elastic  stress- 
strain  relations  apply. 
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c)  In  the  neutral  state  F  vanishes  as  in  case  a,  but  neither  energy 


dissipation  nor  permanent  deformation  occurs  and  the  elastic 
stress-strain  relations  apply.  In  the  present  problem  neutral 
regions  will  be  encountered  in  wnich  neither  the  stress  nor  the 


strain  changes, 


0. 


For  the  purpose  of  this  paper  it  is  convenient  to  combine  elastic  and  neutral 
regions,  jointly  to  be  called  "nondlssipative",  as  opposed  to  plastic  regions, 
where  X  >  0,  indicating  that  energy  is  dissipated.  In  the  nondlssipative  regions 
the  changes  in  stress  and  strain  are  governed  by  the  elastic  relations,  while  in 
plastic  regions,  Eqs.  (7)  and  (8)  apply.  Formally,  the  equations  in  nondissipative 
regions  can  therefore  be  obtained  by  substitution  of  X  =  0  into  the  differential 
equations  derived  below  for  the  plastic  regions  and  by  replacing  the  conditions 
F  »  F  »  0  by  the  inequality  (3). 


J.  is  the  first  Invariant  of  stress,  s  ,  s  and  v  ,  v  are,  respectively,  the 
X  x  y  x  y 

stress  deviators,  and  particle  velocity  components  in  the  x-  and  y-directions. 


1 
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Further,  there  are  two  equations  of  motion 


S*x  ,  1  M1  *  nK 

sr  +  3  sr  +  s?‘ p  sr 


^■  +  i  sr +  l£’p  ss* 


(11-12) 


Equations  (11)  and  (12)  and  the  respective  requirements  on  F  and  X  complete  the 
formulation  except  for  initial  and  boundary  conditions. 


In  plastic  regions  additional  equation  F  *  0  applies,  so  that  there  are  a 

total  of  seven  relations  for  the  seven  unknown  quantities  s  ,  s  ,  T,  J..  ,  v  , 

x  y  x  x 

v  and  X  >  0.  In  nondissipative  regions  Eqs.  (11),  (12)  apply,  but  the  function 
X  vanishes  identically,  X  ■  0,  while  F  must  satisfy  either  F  <  0,  or  the  two 
conditions  F  =  0,  F  <  0  simultaneously.  In  the  nondissipative  case  there  are 
only  six  differential  equations  and  six  unknown  quantities.  The  complete  solution 
of  the  problem  is  to  be  obtained  from  the  six  differential  equations  (11),  (12) 
and  the  applicable  relations  on  X  and  F,  subject  to  appropriate  boundary  or 
initial  conditions  at  the  surface  and  at  the  junctions  of  the  as  yet  unknown 
regions: 


1.  On  the  surface,  y  »  0,  a  step  pressure  p  ■  PQH(Vt  -  x)  normal  to  the 
surface  is  applied,  so  that, 


a  * 


-  P, 


(for  Vt  >  X) 
(for  Vt  <  X) 


(11-13) 


while 


T  »»  0 


0 


(n-i4) 


r 


It  is  known  from  a  general  study  of  elastic -plastic  wave  propagation  [8] 
that  the  largest  characteristic  velocity  possible  is  cp  ,  the  velocity 
of  elastic  P-vaves.  All  stresses  and  velocities  must  therefore  vanish 
outside  the  wedge  formed  by  the  loaded  portion  of  the  surface  and  the 
P-front,  which  is  inclined  at  the  angle 


cp_  «  TT  -  sin 


with  the  x-axis,  Fig. 


-1  1 


2G(l»vj 

pTi^ 


(11-15) 


Because  the  steady-state  problem  is  considered,  stresses  and  velocities  can 
net  be  functions  of  x  and  t  separately,  but  must  be  of  the  form  f(x  -  Vt).  For 
the  step  load  p  *  poH(Vt  -  x)  dimensional  considerations  discussed  in  detail  in 
[5]  require  further  that  the  stresses  and  velocities  do  not  depend  on  x  -  Vt  and  y 
separately,  but  must  be  solely  functions  of  the  combination 


x  -  Vt 


(II-16) 


Instead  of  using  5  as  variable,  the  equivalent  but  more  convenient  one 


-  cot"  ’  5 


(11-17) 


is  introduced.  The  angle  cp  is  shown  in  Fig.  3' 


Noting 


- 5—  ,  one  obtains  the  relations 

sin  $ 


9  _  Id  _  sin  <p  d 
<§x  ~  yd?--  y  dep 


(II-18) 


9  ?  d  sin  2?  d 

57*  -  y  Mzx~~%r  * 


(II -19) 


9  V  d  V  .  2  d 

St  "  -  y  dS  “  y  sln  *  S? 


(11-20) 


Ml-* 


i 


Substitution  into  Eqs.  (11),  (12)  reduces  these  equations  to  a  set  of  ordinary 
simultaneous  differential  equations  in  the  single  independent  variable  ep. 

Because  of  the  manner  of  solution  to  be  employed  the  unknowns  s  ,  s  and  T 

x  y 

are  replaced  by  three  other  dependent  variables,  s^  ,  Sg  and  9, 


2  2 

s  *  s..  cos  8  +  s_  sin  8 

X  -L  d 

2  2 

s  =  sp  cos*-  0  +  s1  sin  6 
yd  x 

T  =  (s.  -  s_ )  sin  9  cos  8 

JL  d 


(11-21) 

(11-22) 


(11-23) 


s^  and  Sg  are  the  two  principal  stress  deviators,  while  8  is  in  the  angle  between 
the  direction  of  and  the  horizontal,  Fig.  3*  Introducing  further  the  angle  y 
between  the  direction  of  and  the  position  vector.  Fig.  3, 


y  *  <p  -  6 


(11-24) 


the  six  differential  equations,  (11),  (12)  become  finally 


•  2  n  '  2  1 

s1  cos  8  +  b2  sin  8  -  (s1  -  8g)  8  sin  28  + 

l_2v  *  2  2 

+  J(l-fv)  J1  +  L^sl  cos  9  +  s2  sin  9) 


V  Yx 


(11-25) 


s^  sin 


i2  8  +  Sg  cos2  8  +  (s^  -  Sg)  8  sin  20  + 

l-2v  „•  „ ,  ,  2  «  2«, 

3(l+vJ  J1  +  L^sl  6in  9  +  s2  cos  9)  = 


2G  ' 

Y  Vy  C0t 


(11-26) 


1-2V  •  2G  '  .  '  . 

I+v”  Jl  =  T  ^Vy  COt  *  -  V 


(11-27 ) 


(81  "  S2^  Sin  20  +  2^S1  "  s2^  ®  cos  29  +  L^si  "  S2 )  sin  29  - 


G  ,  '  .  » 

y  (vx  cot  <P  -  v  ) 


(11-28) 


f 


_  L 

-  9  - 

>  t  i 

sx  cos  0  sin  v  +  sg  sin  0  cos  y  -  (s^  -  sg)  8  cos  (y  -  0)  + 

1  '  • 

+  ^  Jx  sin  9  =*  -  pV  v^  sin  «p  (11-29) 

i  t  i 

s^  sin  0  sin  y  -  sg  cos  0  cos  y  +  (s^  -  3g)  0  sin  (y  -  0)  - 

1  1  * 

-  J.  cos  cp  =  -  pV  v  sin  9  (11-30 ) 

o  y 

Primes  indicate  differentiation  with  respect  to  cp  and  the  function  L  is  related 
to  X  j 

L  *  — ^p—  X  (11-31) 

V  sin  cp 

The  function  L  is  subject  to  the  same  conditions  as  X,  i.e.  L  >  0  in  plastic 
regions,  L  =  0  elsewhere. 


The  expression  for  the  plastic  potential  in  these  variables  is 


2  2  2 
p  *  si  +  S1  *2  +  s2  -  k 


(II -32) 


In  plastic  regions  Eq.  (6)  requires  F  =  0,  giving  the  additional  differential 
equation 


(2bx  +  sg)  sx  +  (sx  +  2sg)  sg  ~  0 


(11-33) 


t  t 

The  unknowns  v  and  v  can  be  eliminated  from  Eqs.  (25)-(30)  without 
x  y 

differentiation.  Using  the  symbol 

OV2  p 

X  «  X(cp)  =  sin  (p 


(H-34) 


the  operations 


(11-35) 


Eq.  (40a)  «  Eq.  (25)  +  Eq.  (26)  -  Eq.  (27) 

Eq.  (40b)  *  -  X  Eq.  (27)  +  sin  cp  Eq.  (29)  -  cos  9  Eq.  (30) 


(11-36) 
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Eq.  (40c)  =  X  {sin  28  [Eq.  (26)  -  Eq.  (25)]  +  cos  20  Eq.  (28)}  + 

+  cos  (y  -  8)  Eq.  (29)  -  sin  (Y  -  ®)  Eq*  (30)  (11-37) 

Eq.  (40d)  =  X  {cos  2©  [Eq.  (26)  -  Eq .  (25)]  -  sin  20  Eq.  (28)}  + 

+  sin  (y  -  8)  Eq.  (29)  -  cos  (y  -  0)  Eq.  (30)  (II-38) 

Eq.  (UOe)  *  Kq.  (33)  (II -39) 


lead  to  the  following  set  of  five  differential  equations  valid  in  plastic  regions: 


— 

" 11,1 

1 

1 

-O 

0 

(S1  +  S2J 

1 

S1 

2 

sin  Y 

2 

cos  y 

1-3X  ffi) 

-  sin  2y 

0 

1 

S2 

\  sin  2y 

\  sin  2y 

sin  2y 

2X  -  1 

0 

1  r' 

3  J1 

sin2  y  -  X 

2 

X  -  cos  y 

-  cos  2y 

0 

-  X(s1  -  s2) 

(S1  -  s2)0 

2sl  +  S2 

S1  +  2s2 

0 

0 

0 

L 

(II-UO) 


The  system  of  equations  for  nondissipative  regions  consists  of  the  first 
four  Eqs.  (40)  without  the  terras  containing  L,  i.e. 


I 


-  11  - 


-  1 

-  i 

1-2V 

1+v 

2 

sin  v 

2 

cos  v 

1-2V 

1-3X 

sin  2y 

sin  2y 

2  sin  2y 

sin2  v  -  X 

v  2 

X  -  cos  y 

-  cos  2y 

0 

t 

S1 

-  sin  2y 

t 

S2 

-  2(1-2X) 

1  t' 

3  J1 

0 

(Sx  -  S2 

mrmm 

MM 

(II-Ul) 
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III  SOLUTIONS  FOR  INDIVIDUAL  REGIONS. 

As  a  first  step  towards  the  construction  of  overall  solutions,  expressions 
for  individual  regions  must  be  derived.  The  latter  will  be  combined  in  Section  IV 
to  find  the  solution  for  the  entire  domain. 


a )  Nondisaipative  Regions. 


Equations  (II-41)  are  linear  and  homogeneous  so  that  the  derivativea 

f  f  f  » 

of  the  stresses  ,  sp  ,  and  the  value  (e^  -  ep)  9  vanish,  unless 
the  coefficient  matrix  in  Eqs.  (11-41)  is  singular,  requiring 

X(1-2X)  (1  ♦  (l-2X)(l-2v)]  .  0  (III-l) 


Equation  (1)  has  two  significant  roots, 


and 


*P 


1  -v 
l-2v 


1 

2 


(III-2) 


(HI-3) 


Substitution  of  the  two  roots  Xp  and  Xg  into  (11-34)  furnishes  the 
two  locations 


where  cp 


-1  cp 
n  -  sin  1  f 


tt  -  sin 


-1  CS 


(III-4) 

(HI -5 ) 


rS  V 

Cg  are  the  velocities  of  P-  and  S-waves,  respectively.  In 


all  locations  tp  «>  cpp  or  cpg  the  values  ,  (e^  -  «g)  9  vanish 

so  that  in  nondissipative  regions  the  stresses  must  remain  constant 
except  at  the  locations  tpp  and  «pg  .  The  latter  being  the  potential 
locations  of  elastic  P-  and  S-shock  fronts,  respectively,  it  is  known 
that  discontinuities  in  stresses  and  velocities  may  occur  at  these 
locations  and  may,  therefore,  be  part  of  the  complete  solutions  to  be 
constructed.  The  following  pertinent  details  will  be  required  subse- 


r 


quently . 
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(1)  The  P-front. 


Designating  the  discontinuous  changes  in  the  various  quantities 
at  the  front  by  the  symbol  A,  the  discontinuities  in  the  stresses 

u  az  (normftl  and  tangential  to  the  front,  respective.'' y)  and 
in  the  component  of  the  velocity  (normal  to  the  front)  are  proportional 


to  . 


(III-6) 


No  other  discontinuities  can  occur  in  this  location. 


The  changes  bo^  and  bo ^  are  of  course  limited  by  the  yield  condition 
F  <  0  which  must  be  satisfied  on  either  side  of  the  front. 


In  the  actual  solution  a  P-front  will  be  encountered  only  when  the 
region  ahead  of  the  front  is  stressless  and  at  rest.  The  normal  to  the 
front  is  then  a  principal  direction  for  the  stresses  behind  the  front, 
so  that  v  *  0  or  ^  .  Selecting  ,  corresponds  to 


Y 


(III -7) 


The  corresponding  values  of  the  other  quantities  of  interest  behind 


the  shock  front  are 

2(1 

e  ’  3  ’  '1  ■  3(1 ’  “2 


(1-2V)0  (1+Y)© 

ITITvr  '  J1  ’  “T^T“  (HI -8) 


subject  to  the  limitation 


l-2v 


( III -9 ) 


imposed  by  the  yield  condition. 
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(2)  The  S -front. 


At  an  S-front  discontinuities  occur  only  in  the  shear  stress 
t^  ■  T^  a  t  and  in  the  tangential  velocity  .  The  change  in  velocity 

is  proportional  to  At 


(III-10) 


In  addition,  the  yield  condition  F  <  0  must  again  be  satisfied  ahead 
of  and  behind  the  front. 


The  relations  between  the  state  of  stress  on  either  side  of  an 
S-front  in  terms  of  At  and  of  the  variables  s^  ,  s^  ,  B  and  y  can  be 
obtained  in  such  routine  manner  that  only  one  detail  used  in  Section  IV 
is  presented  here. 


It  is  possible  for  an  S-front  to  occur  between  two  neutral 
regions,  i.e.  regions  of  constant  stress  for  both  of  which  the  yield 
condition,  F  ■  0,  is  satisfied.  In  this  special  case  the  quantities 
3,  ,  a  and  8g  have  no  discontinuity  at  the  front,  only  the  direction 

of  the  principal  stress  changes.  The  values  of  the  angles  y,  V  ahead 
of  and  behind  the  front,  respectively,  are  complementary 

V  ■  n  -  y  (III-ll) 

as  shown  in  Fig.  (5 ). 


b  )  Plastic  Regions. 

In  such  regions  Eqs.  (11-40)  apply.  They  are  linear  and  homogeneous 

I  f 

in  the  value3  ,  Sg  ,  etc.,  anc.  may  be  satisfied  by 

Ilf  f 

S1  “  S2  ”  J1  “  (S1  "  s2^  9  “  L  =  0 


( III-12 ) 
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I 

f 

I 

I 

I 

I 

I 

i 

I 

l 

I 


However,  L  -  0  implies  X  -  0  which  violates  Eq.  (II-9).  It  follows 
that  in  plastic  regions  the  determinant  of  Eqs.  (II-UO)  must  vanish,  giving 
the  determinantal  equation 


(®!  +  s2)2  (bg  +  bxb3) 


0 


where 


b1  -  2[1  +  (l-2v)(l-2X )J 
b2  -  0  cos  2y  +  (l-2X)(l-2v) 
b3  -  (1+V)(1-2X)  -  02X 


and 


(Hi-13 ) 

(III-lU) 
(III-15 ) 
(III-16) 

(IH-17) 


If  this  value  is  substituted  it  will  be  seen  that  Eq.  (13)  is  a  homogeneous 
quadratic  expression  in  s^^  . 


Due  to  the  vanishing  of  the  determinant  only  four  of  the  five 
Eqs.  (II-40)  are  Independent .  By  definition  L  must  not  vanish,  so  that 

r  *  r  r 

s i  >  °2  >  8111(1  ®  can  always  be  expressed  in  tenns  of  L, 

,  (3-0)  bu  (C]L  +  bs)  L 


3b, 


(3+0)  b4  (bx  +  s2)  L 


3b, 


,  sin  2y  b,L 

9  3  5Ti72x7  bl 


Ji  ■  l>2  -  I  * b  J  — 


(s,  +  E  )  L 


(III-lfl) 

(III-19) 

(III-20) 

(III -21 ) 


Velocities  and  accelerations  may  be  obtained  from  the  relations 


where 


V  sin  cp  (s.  +  s  )  L 


v  = 
x 


v  - 

y 


[b2S  sin  (2y  -  <P)  -  ?o3  sin  <p]  (m-22) 


-  V  sin  tp  (s1  +  s2)  L 


2G(l-2X  )Tb~ 


|^b2@  cos  (2y  -  9)  +  2b  ^  cos  9]  (III-23) 


vx  =  1  sin  9  v’ 


(III -24 ) 


v  = 


V  1 

—  sin  cp  v 

R  y 


(III -25) 


b,  «  (1+v)  con  2y  +  gX(l-2v) 


(111-26) 


Since  Eq.  (13)  must  remain  valid  throughout  a  plastic  region,  it  may 
be  differentiated  with  respect  to  9.  This  leads  to  an  expression  which 
contains  the  first  derivatives  of  the  stresses  linearly,  so  that  substitution 
of  Eqs.  (l8)-(20)  furnishes  a  linear  equation  for  L.  Its  solution  gives  L 
as  a  function  of  3,  Y  and  of  the  position  angle  9, 


3bp(l-2X) 
— £ - -  { 

4  sin2  9 


X  sin  29  [4(l-2V)(bg  +  b^)  +  b]L(32+  2  +  2v)]  +  4b^  sin  2y  sin2  9 
(3+32)(l-2X)  b4  (b2  cos  2y  -  b^X)  +  3bgb3  sin2  2y 


(HI-27) 

rtt  t 

The  values  of  the  derivatives  s^  ,  Sg  ,  and  6  can  be  obtained  by 
substitution  of  Eq.  (27)  into  Eqs.  (l8)-(25)* 


In  principle  Eqs.  (l8)-(25)  permit  the  numerical  determination  of  the 
values  of  stresses  and  velocities  in  the  interior  of  a  plastic  region  by 
quadratures  if  the  values  on  one  boundary  of  this  region  are  known.  The 
starting  values  must  inherently  satisfy  the  yield  condition,  F  =  0,  and 
the  determinantal  equation  (13 ) *  Further,  throughout  the  plastic  region, 


(ni-28) 


-  IT  - 


L  >  0 


must  be  satisfied. 


Discontinuities  (Plastic  Shock  Fronts). 


It  is  known  that  in  transient  problems  one,  but  just  one  type  of  plastic 
shock  front  can  propagate  in  the  elastic-plastic  material  considered  here  [8j. 
However,  such  a  front  can  exist  only  in  locations  where  the  normal  to  the 
front  l<es  in  the  direction  of  one  of  the  principal  stresses,  while  the 
other  two  are  equal  and  where  the  yield  condition  is  satisfied.  The  velocity 
of  propagation  of  the  front  is 


(III-29) 


where  K  "  G  is  the  bulk  modulus.  The  discontinuity  is  restricted 

to  the  particle  velocity  v^  normal  to  the  front  and  to  the  first  invariant 


.  The  change  must  have  the  same  sign  as 


AJ, 


>  0 


(111-30) 


The  other  conditions  stated  define  the  values  of  y  and  P 


TT 

Y=  2 


(III-31) 


A  discontinuity  traveling  in  real  space  with  the  velocity  c,  Eq.  (29), 
can  occur  in  the  steady-state  problem  only  in  the  location 


9  *  tt  -  sin""1  (~) 


(III- 32) 


The  corresponding  value  of  X  is 


X  = 


1+V 

3U-2V) 


(ni-33) 


The  denominator  in  Eq.  (27)  vanishes,  as  ejected,  for  these  values  of 
y,  P  and  X. 


1 
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The  possibility  of  the  occurrence  of  this  discontinuity  (plastic  shock 

front)  must  be  considered  when  constructing  the  complete  solutions  in 

Section  IV.  It  was  actually  found  that  no  such  shocks  occur,  except  in  the 

p0 

limit,  V  -*■  «  .  However,  for  large  values  of  the  parameter  defining  the 
surface  load,  the  solutions  come  extremely  close  to  the  singular  values 
representing  a  shock,  so  that  computing  difficulties  occur. 

d)  Asymptotic  Solutions  near  Singularities. 

As  stated  in  the  previous  paragraph  numerical  difficulties  in  the 
vicinity  of  9  »  $  will  make  the  procedure  for  integration  of  Eqs.  (18-27) 
outlined  in  subsection  b  unsuitable  if  the  values  of  3  and  y  are  sufficiently 
close  to  those  for  a  plastic  front,  3  »  3>  Y  *  • 

To  establish  the  behavior  of  the  solution  of  Eqs.  (18-27)  in  such  cases, 
let 


3  *  3  +  A  (111-34) 

9  «  $  +  t 

where  1),  A  and  «  are  small  quantities.  Substitution  of  Eqs.  (34)  into  the 

t  f 

deterrainantal  equation  (13  )^  Eq.  (20)  for  9,  and  an  equation  for  ^  *  A 
obtained  from  a  combination  of  Eqs.  (17 ),  (18)  and  (.19),  provides  three 
simultaneous  first  order  differential  equations  for  the  three  unknowns 
1),  A  and  L,  subject  to  the  inequality  L  >  0.  These  three  nonlinear 
differential  equations  are  of  course  equivalent  to  the  equations  in  sub¬ 
section  b  and  no  easier  to  treat.  However,  A,  7)  and  c  being  small  quantities, 
approximate  equations  may  be  obtained  by  neglecting  higher  order  terms. 
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Various  degrees  of  approximation  are  possible  depending  upon  tne  number  of 

higher  order  terras  of  1),  A  and  t  which  are  retained.  Appendix  A  presents 

the  resulting  set  of  equations  when  the  leading  and  the  next  terms  in  each 

variable  are  retained.  In  [5j  the  comparable  set  of  equations  is  given  when 

only  leading  terms  of  each  unknown  are  retained.  Together  vith  the  appropriate 

bJ1 

expression  for  the  approximate  equations  may  be  used  when  integrating  in 
the  vicinity  of  9  *  9. 

p0 

To  obtain  results  for  high  values  of  -j-  the  integration  of  the  original 
differential  equations  (18-27)  was  stopped  when  accuracy  troubles  developed 
and  the  integration  was  continued  using  the  approximate  equations  derived  in 
Appendix  A.  This  process  was  successful  for  all  numerical  cases  considered. 

The  results  obtained  in  this  manner  in  the  vicinity  of  the  singular  point  approach 
simple  asymptotic  expressions  derived  in  [5J  for  v  yf  g-  .  This  is  shown  for  typical 
cases  in  Figs.  15,  17  and  l8  of  this  reference  and  provides  a  check  on  the  results 
of  the  complicated  analysis  in  Appendix  A. 

For  the  special  case  v  =*  g  an  asymptotic  expression  valid  for  9  <  9g  ■  9  is 

given  in  Appendix  B.  As  demonstrated  on  a  typical  case  in  the  same  Appendix,  it 

1  Po 

is  possible  to  construct  solutions  for  v  =  g  for  high  values  of  •£-  by  using  this 
asymptotic  solution  for  9  <  9-  and  a  numerical  integration  of  the  original  differ- 

O 

ential  equations  (18-27 )  with  double  precision  for  9  >  9g  . 


IV  CONSTRUCTION  OF  SODJTIONS 


In  Section  IV  a  number  of  partial  solutions  were  obtained  from  which  the 
solution  of  the  complete  boundary  value  problem  is  now  to  be  constructed. 

Section  Ill-b  gives  tne  differential  equations  for  the  determination  of  the 
stresses  and  velocities  in  plastic  regions;  Section  Ill-a  indicates  that  all 
unknowns  in  nondissipative  regions  are  constants,  except  for  discontinuities  of 
a  prescribed  nature  at  the  locations  9g  and  9p  .  In  addition,  as  discussed  in 
Section  III-c,  a  shock  front  with  plastic  deformation  may  occur  at  the  location  9. 


In  Section  II,  boundary  conditions  and  additional  requirements,  which  the 
solution  must  satisfy,  were  formulated  and  discussed.  Equations  (11-13  and  lb) 
for  the  prescribed  surface  load  in  terms  of  the  variables  s^  ,  J.,  ,  y  and  9  require 
either 


Jl(") 

81(TT)  +  - - -  *  -  PQ 


(IV-la) 


or 

JX(TT) 

a2(T?)  +  *  -  pQ  ,  »(n)  -  0 


(IV-lb ) 


A  further  boundary  condition  requires  that  all  quantities  must  vanish  for 
9  <  9p  ,  Eq.  (11-15 )•  This  condition,  in  conjunction  with  the  fact  that  a 
plastic  region  or  a  plastic  shock  can  exist  only  in  locations  where  the  yield 
condition  is  satisfied,  permits  the  conclusion  th8,t  the  change  in  stress  from 
vanishing  values  for  9  <  9p  to  nonvanishing  values  must  be  nondissipative. 
However,  in  nondissipative  regions  the  stresses  are  constant,  except  for  dis¬ 
continuities  at  9  «  9 p  or  9  »  9g  .  A  solution  in  which  plastic  deformations 
occur  at  all  can  therefore  start  only  in  one  of  the  two  ways  described  below. 
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Case  1.  Discontinuities  occur  at  the  P-  and  S«>fronts,  where  the  discontinuity 
at  cp_  satisfies  the  inequality 

ST 


a  L(+)l  < 

1  LfP  J  1-2  V 


(IV-2) 


while  the  discontinuity  AT  at  cpg  is  of  such  magnitude  that  the  yield  condition 

F  [<4+)]  -  0  (IV-3) 

is  satisfied.  The  symbols  ^  or  ^  indicate  approach  from  above  or  below, 
respectively. 


Case  2.  A  discontinuity  occurs  at  the  P-front,  described  by  Eqs.  (III-7>8), 

(  +  ) 

so  that  the  yield  condition  is  satisfied  for  qp  *  cpp  ,  i.e. 


(IV-4) 


In  Case  1,  plasticity  may  occur  only  in  locations  >  epg  ,  while, 
in  Case  2  it  may  occur  already  for  cp,0  <  9  <  «pg  . 


As  a  next  step  in  the  search  for  solutions  it  is  helpful  to  consider  the 

Pc 

latter  as  function  of  the  nondimensicnal  surface  pressure  -j-  ,  while  Poisson's 

ratio  V  and  the  velocity  V  are  considered  constant.  For  sufficiently  small 

Po  Po 

values  of  —  the  solution  must  be  entirely  elastic,  but  as  increases  plastic 

regions  must  occur  and  should  form  a  gradually  changing  pattern.  Based  on  [1] 

P»T 

one  can  find  the  limiting  value  ~  ,  above  which  entirely  elastic  solutions  are 
no  longer  possible. 


The  elastic  solution,  shown  in  Fig.  6,  has  two  discontinuities  at  «p?  and  «pg 

with  regions  of  constant  stress  between  and  sad  between  «pg  and  the  loaded 

y 

surface  tp  »  17*  Depending  on  the  values  of  the  parameters  v  and  — -  ,  the  yield 

P 

condition  may  be  reached  in  either  of  the  two  regions  resulting  in  different 


1 


\ 


-  22  - 


expression  for  y  .  If  the  region  9  >  <pg  controls ; 


E 


[ 


3N2  -  3N  COB  2q>g  +  (1-V+^)  COB2  2?s 


'] 


(IV-5a) 


where 


N  -  |  [cob  2sps  +  (1-2V)  cos  2(epg  -  <pp)] 


(IV~5b) 


while 


[ 


(1-2 v)2  COB2  2cp, 


-1 


(lv-6) 


if  the  region  cp  <  cp  controls.  The  decision  which  region  controls  can  be  made 

o 

by  comparing  the  values  given  by  Eqs.  C>a)  and  (6),  the  smaller  one  controlling. 

y 

Designating  as  Range  I  the  combination  of  values  V  and  — -  where  Eq.  (5a)  controls, 

cp 

one  finds  that  in  this  range 

(0  >  (w"7) 


The  remainder  of  the  range  will  be  designated  as  Range  II.  Both  ranges  are  shown 

Pr 

in  Fig.  7*  The  limiting  values  y  are  shown  in  Fig.  8  for  several  values  of  V  as 
V 

function  of  —  . 

CP 


PE 

If  the  surface  load  exceeds  the  value  y  by  a  sufficiently  small  amount  the 
elastic-plastic  solution  should  differ  only  slightly  from  the  elastic  one,  which 
can  be  used  to  predict  the  character  of  the  solution.  Because  the  situations 
differ,  the  Ranges  I  and  II  must  be  discussed  separately. 


a)  Solutions  in  Range  I. 

PE 

In  this  case  y  is  given  by  Eq.  (5a)  and  the  yield  stress  in  the  elastic 
solution  is  reached  only  in  the  region  cp  >  «pQ  *  The  discontinuity  satisfies 


L 
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then  the  inequality  (2)  and  continuity  requires  that  this  inequality  will  still 


PE 


apply  for  a  range  of  values  >  -g*  >  ~  >  where  is  a  limiting  value,  not  yet 
known.  In  this  range  the  start  of  the  solutions  will  be  according  to  Case  1  and 
plasticity  can  therefore  occur  only  in  the  region  qp  >  «pg  . 


Using  an  indirect  approach,  the  determination  of  plastic  solutions  for  values 
Pn  PT 

-1-  <  -~  begins  with  the  selection  of  e.  pair  of  starting  values  a.  and  At  near  those 

K  K  «*• 

for  the  limiting  elastic  case.  Experience  and  continuity  considerations  indicate 

that  this  value  should  be  larger  than  ©^  corresponding  to  pg  given  by  Eq.  (5e). 

A  plastic  region  can  start  only  at  a  point  which  is  a  root  of  Eq.  (III-13). 

Inequalities  derived  in  [5]  indicate  that,  for  a  given  state  of  stress,  this 
equation  has  only  two  roots  in  the  meaningful  range  of  <p,  namely  "■  <  cp  <  tt,  subject 
to  the  following  bounds: 


<Pp  <  <PA  <  <PS  >  <Pg  <  *b  <  "  (IV-8) 

Starting  integration  at  ep^  «  cp^  the  solution  in  the  interior  of  the  plastic  region 
is  determined  from  Eqs.  ,III-l8  to  27).  The  plastic  region  can  be  extended  as 
long  as  Eq.  (III-27)  gives  values  L  >  0,  but  the  plastic  region  may  be  terminated 
at  will  at  any  earlier  location  <pp  .  The  solution  for  cp  >  <p2  is  then  nondissipative, 
i.e.  all  quantities  are  constant.  If,  therefore,  in  the  process  of  forward  in¬ 
tegration  a  value  ®  is  encountered  which  satisfies  Eq».  (1),  the  plastic  region 

P0 

is  terminated  and  a  solution  for  one  value  of  the  surf  ace  pressure  -r-  has  been 
obtained.  Repeating  this  process  with  gradually  increasing  starting  values 
the  whole  spectrum  of  values  satisfying  the  inequality  (2)  can  be  explored. 

Solutions,  if  any,  obtained  in  this  manner  will  have  the  configuration  *hovn  in 
Fig.  9,  i.e.  discontinuities  at  cpp  and  epg  and  a  plastic  region  tt  >  >  «px  >  «pg  . 

There  will  be  an  elastic  region  of  constant  stress  from  to  <pg  and  two  neutral 
regions  to  either  side  of  the  plastic  one. 


1 
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C.  1; 


One  can  proceed  in  a  similar  manner  when  the  solution  begins  at  9  *  <pp  as 

indicated  in  Case  2.  Starting  with  a  value  according  to  Eq.  (3)>  the  yield 

condition  is  satisfied  for  any  value  ep  >  cpp  ,  so  that  the  determinantal  equation 

(III-13)  according  to  Eqs.  (8)  now  has  two  roots,  9^  >  9g  >  v^ich  plastic 

regions  may  start.  Both  roots  must  be  explored.  If  the  larger  one,  9  >  9g  , 

leads  to  a  solution,  it  has  a  configuration  as  shown  in  Fig.  10a.  Starting, 

alternatively,  with  the  smaller  root,  <pp  <  9^  *  9^  <  9g  >  several  possibilities 

are  to  be  investigated.  The  integration  may  be  continued  as  long  as  L  >  0  to  see 

if  a  value  0  «  0  or  |  can  be  reached.  The  configuration  of  such  a  solution,  if 

any,  is  shown  in  Fig.  10b.  Alternatively,  the  plastic  region  can  be  terminated 

at  will  at  a  point  9g  <  «pg  where  0  jl  0,  ~  .  The  plastic  region  will  then  be 

followed  by  a  neutral  one  for  values  9  >  9^  .  The  inequalities  (8)  on  the  roots 

of  Eq.  (III-13)  indicate  that  there  is  just  one  more  root  9^  >  9g  ,  when  a  second 

plastic  region  can  begin.  Starting  integration  at  this  point  may  lead  to  a 

terminal  location  9^  ,  where  0  «  0  or  -  .  The  configuration  of  such  a  solution, 

if  any,  Fig.  10c,  contains  a  P-front  and  two  plastic  regions,  separated  by  three 

neutral  regions.  There  are,  however,  further  possibilities.  The  neutral  region 

9  >  9g  which  follows  the  first  plastic  one  may  be  terminated  at  9g  by  an  elastic 

change  in  shear,  At,  which  is  restricted  in  sign  and  intensity  by  the  yield 

condition.  If  AT  is  such  that  F  |jpg+^J  <  0,  "the  region  9  >  9g  becomes  elastic. 

This  might  permit  values  0  b  0,  jj  at  the  surface,  the  corresponding  solution 

having  the  configuration  of  Fig.  11.  Finally,  the  important  case  must  be 

considered  where  the  value  of  At  is  such  that  F  [Vg+^j  vanishes  again,  a  situation 

discussed  in  Section  III-a-2.  In  the  latter  case  there  is  again  a  neutral  region 

for  9  >  9  ,  which  can  be  followed  by  a  plastic  region  because  Eq,  (III-13)  has 

a  root  90  >  90  giving  a  starting  location.  The  configuration  of  a  solution 
3  5 

obtained  in  this  manner  is  shown  in  Fig.  12. 
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b )  Solutions  in  Ranpje  II. 


PE  P0  pv 

In  this  range  is  given  by  Eq.  (6)  so  that  in  the  limiting  case  y  - 


yield  is  just  reached  in  the  region  qPp  <  cp  <  cpg  .  The  discontinuity  at  the 

P-front  must  therefore  satisfy  Eq.  (4),  which  will  also  hold  for  neighboring 

Pq  pe 

elastic-plastic  solutions  where  —  exceeds  ~  slightly.  These  solutions  will 

therefore  start  at  q>  =*  cp_  according  to  Case  2.  In  the  limiting  solution  for 
P0  PE 

IT  *  IT  reSion  >  ^5  below  yield  and  continuity  requires  this  to  hold  in 

neighboring  elastic-plastic  solutions,  so  that  the  plastic  region  must  lie  in  the 

range  «pp  <  9^  <  <  cpg  •  The  construction  of  solutions  begins  exactly  as  for 

P0  PE 

IT  >  IT  *  For  each  terralnal-  Point  cp2  the  strength  of  the  discontinuity  AT  at  the 
shear  front  is  determined  by  the  requirement  that  8  ■  0  or  |  subject  to  the 
limitation  F 


[#*]  <  0.  When  the  required  value  At  violates  this  condition  a 
second  plastic  region  for  cp  >  <pg  is  needed,  i.e.  the  configurations  shown  in 
Figs.  10c  and  12  are  to  be  investigated. 


c )  Alternative  Solutions  and  Considerations  of  Uniqueness  and  Existence. 

In  the  absence  of  a  uniqueness  theorem  it  is  vital  to  demonstrate  that 
configurations  other  than  those  in  Figs.  9  to  12  can  not  lead  to  solutions. 
According  to  [5],  Eq.  (111-3,3)  has  for  a  given  state  of  stress  one  root  cp,  no 
more  no  less,  in  each  of  the  intervals  <p_  <  cp  <  <p  and  <p  <  cp  <  tt.  If  a  plastic 

a  D  D 

region  ends  at  a  location  in  one  of  these  intervals,  the  state  of  stress  in 
the  remainder  of  the  interv  tl  for  tp  >  cp^  is  necessarily  neutral  and  uniform  and 
equal  to  the  one  at  the  terminal  point  of  the  plastic  region.  Cp^  is  therefore 
the  only  solution  of  Eq.  (III-13)  for  this  state  of  stress  in  the  particular 
interval  and  no  more  than  one  plastic  zone  can  therefore  occur  in  any  interval. 
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In  Section  III~3  the  possibility  of  discontinuous  plastic  shock  fronts  has 
been  indicated  and  their  occurrence  must  be  considered.  It  has  been  shown  in  [5] 
that  for  finite  values  of  V  no  plastic  shock  front  can  occur  and,  that  there  can 
be  no  more  than  one  plastic  region  in  each  of  the  intervals  cpp  <  cp  <  <pg  , 

<pg  <  cp  <  7T.  While  discontinuous  plastic  shock  fronts  can  not  occur,  values  of 
Y  and  3,  where  the  conditions  (III-31)  are  nearly  satisfied,  are  encountered. 

The  asymptotic  behavior  of  the  solution  near  <p  *  $  in  such  cases  was  studied  in 
Section  Ill-d  and  details  are  given  in  Section  5  of  [5].  Therefore,  combined  with 
elastic  discontinuities  at  the  P-  and  S~fronts,  only  the  limited  number  of 
configurations  shown  in  Figs.  9  to  12  are  possible. 

The  numerical  analysis  by  digital  computer  was  set  up  to  investigate  all 
possible  alternatives,  i.e.  the  configuration  according  to  Fig.  9  if  the  starting 
value  <3^  satisfies  Eq.  (2)  and  any  of  the  alternatives  shown  in  Figs.  10  to  12  if 
cr1  satisfies  Eq.  (3).  While  none  of  the  configurations  shown  in  Figs.  lOa-c  ever 
furnished  a  solution,  no  general  proof  permitting  elimination  of  these  cases  is 
available. 


In  Range  I  solutions  which  start  according  to  Case  1,  have  the  configuration 

of  Fig.  9.  For  fixed  values  of  V  and  V,  these  solutions  form  a  family  which 

depends  on  one  parameter,  the  selected  starting  value  a,  >  a  .  It.  was  found 

p  1  IE 

that  the  surface  load  -rr  increases  raonotonically  with  a,  until  the  limit.  Ea.  (4) 

iv  x  .... 

pr 

for  a.  is  reached,  which  leads  to  a  limiting  value  of  the  surface  load  ~  . 

P  k 

However,  no  analytical  proof  of  the  monotonic  increase  of  is  available. 


The  solutions  found  for  Range  I,  which  start  according  to  Case  2,  had  always 
the  configuration  shown  in  Fig.  12.  These  solutions  also  depend  on  one  parameter, 
viz.  the  stopping  point  of  the  plastic  region  between  cp_  and  cp_  .  If  tf.  is 
selected  only  slightly  larger  than  cp^  ,  the  solution  must  obviously  be  very  close 
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Po  ~  PL 

to  the  limiting  one  for  Case  1,  so  that  in  such  a  case  there  is  a 

smooth  transition  from  the  configuration  according  to  Fig.  9  to  that  of  Fig.  12. 

P0 

The  numerical  analysis  indicated  that  the  surface  load  -j-  increases  monotonically 
with  <pp  .  As  cp  approaches  a  limiting  value  the  surface  load  goes  to  the  limit 

p 

ro 

-r-  -*■  «  ,  for  reasons  explained  in  Section  Ill-d. 


In  Range  II,  only  solutions  which  start  according  to  Case  2  were  found,  their 

configurations  being  as  6hown  in  Figs.  11  and  12.  Figure  11  applied  for  values 
pE  PQ  PL  PL 

T  <  TT  -  T  w^ere  ~  is  a  bound.  The  corresponding  family  of  solutions  depends  on 

pL 

the  stopping  point  9^  of  the  plastic  region.  The  bound  is  reached  when  the 

P 

elastic  region  for  cp  >  cpg  becomes  neutral.  For  larger  values  of  ~  Fig.  12 
applies  and  all  statements  made  in  Range  1  for  this  case  apply. 


In  Range  I  as  well  as  in  Range  II,  combination  of  all  solutions  obtained 

p0  pE 

numerically  furnished  one,  and  only  one  solution  for  each  value  of  Y*  >  T  * 
However,  no  general  proof  is  available  that  this  must  be  so.  Existence  and 
uniqueness  of  the  solutions  obtained  must  therefore  be  demonstrated  for  each 

V 

combination  of  values  V  and  —  by  actual  computation  of  the  families  of  solution 

CP 

according  to  the  configurations  shown  in  Figs.  9  to  12. 


V  RESULTS  AND  CONCLUSIONS. 


a)  Numerical  Results  for  the  Stresses. 

For  the  numerical  integration  of  the  simultaneous  differential  equations 

(III-18  to  21)  in  plastic  regions,  a  Runge-Kutta  forward  integration  scheme  of 

Fourth  order,  [7],  was  used.  Computations  were  programmed  in  FORTRAN  for  an 

IBM  7090  digital  computer.  Results  for  the  stresses  are  given  in  Tables  1  to  l6, 

(to  be  used  in  conjunction  with  Fig.  13),  for  all  combinations  of  the  parameters 

v  =  0,  0.125,  0.25,  0.35  and  —  =  1.25,  1.5,  2.0,  4.0,  for  five  different  values 
P«  PG  Cp 

.  The  values  -£•  have  been  selected  such  that  there  are  a  sufficient  number 


of 


of  results  to  permit  interpolation  in  each  of  the  configurations  applicable  for 

V 

each  of  the  combinations  of  v  and  —  . 

°p 

p0 

In  general  there  are  three  distinctive  values  of  for  each  case,  which  are 


P™  pt  PH 
called  and 

\  PE  P0  PL 

entirely  elastic.  The  distinctive  value,  >  defines  the  range  Y<  F-T  where 


The  value  is  the  one  up  to  which  the  solution  is 

Pt  ~ 


a  configuration  with  one  plastic  region  applies  according  to  either  Fig.  9  or 
PG  PL 

Fig.  11.  For  -£->“£*  aPPHcai:>le  configuration,  Fig.  12,  contains  two  plastic 
regions. 


As  mentioned  in  Section  III  computational  difficulties  arise  in  the  case 
P0  pL  Pc  p0 

*k*  >  IT  wiien  IT  ^ecomes  l^Ssr  than  about  5*  In  the  integration  the  value  of 

obtained  dependt?  on  the  selected  end  point  op^  of  the  lower  plastic  region-  The 

details  of  the  difficulties  depend  on  the  value  of  Poisson's  ratio,  g  .  For 

v>  |  ,  cpg  remains  smaller  than  9,  but  the  computation  becomes  sensitive  when  9^ 

approaches  9.  For  V  <  g  the  point  9  is  situated  in  the  interior  of  the  upper 

plastic  region,  but  the  computation  in  the  vicinity  of  9  becomes  very  sensitive  to 

small  changes  in  9g  .  For  v  »  g  ,  when  9  ■  9g  ,  both  plastic  regions  approach  9, 

and  the  computation  again  becomes  very  sensitive  to  small  changes  in  90  . 


In  either  of  these  cases  the  sensitivity  is  due  to  a  very  rapid  change,  in  the 


proximity  of  in  the  quantity  .  This  change  is  extremely  rapid,  nearly  a 
shock  front,  while  other  quantities  near  9  approach  limiting  values  smoothly. 


For  values  -g-  where  computational  difficulties  arise,  the  approximate 

differential  equations  (A-5-7)  derived  in  Appendix  A  must  be  used.  One  result, 
pu 

designated  by  —  ,  requiring  this  type  of  analysis  is  given  in  the  tables  for 

V 

each  combination  of  v  and  —  .  The  stresses  ,  <J0  and  J.  for  any  higher  value 

CP  Po  (pH  +  S)  H  H 

of  the  surface  pressure,  say  —  =  -  ,  are  equal  to  +  6  ,  +  6  , 

+  3  ,  respectively,  at  locations  9  >  9,  where  the  superscript  H  indicates 

the  respective  values  for  p  given  in  the  tables.  The  stresses  for  location 

n 

9  <  9,  and  the  quantities  0  and  9  everywhere,  are  approximately  equal  to  those 


for  p„  given  in  the  tables. 

rl 


The  numerical  results  for  the  stresses  are  recorded  in  Tables  l-l6.  The 

notation  "Not  Applicable"  in  the  tables  which  may  appear  in  Regions  B  or  C, 

indicates  that  the  respective  region  does  not  occur  at  all  for  this  value  of  pQ  , 

i.e.  region  A  or  D,  respectively,  or  both  extend  as  far  as  the  S-front.  For 
Po  PH 

values  <  ■—  results  may  be  obtained  by  interpolation  (sometimes  nonlinear), 


po  PH 

while  for  -r-  >  —  the  procedure  stated  in  the  previous  paragraph  is  to  be  used. 

K  K 

At  points  of  transition  in  the  configuration,  i.e.  for  the  loads  pE  and  p^  ,  it 
is  seen  that  the  upper  and  lower  limits  of  one  of  the  plastic  regions  are  equal, 
9n  *  90  ,  or  9,  ■  9]a  indicating  that  there  is  actually  no  plastic  region.  These 
values  are  given  in  the  tables  to  permit  interpolation  for  pressures  exceeding 
the  respective  value  p^,  or  . 


b)  Simplified  Determination  of  Velocities  and  Accelerations. 

The  basic  relations  in  Section  II  permit  the  numerical  determination  of 
stresses  and  velocities  or  accelerations.  The  integration  for  the  stresses  must 
be  actually  carried  out  to  satisfy  the  boundary  condition  on  the  surface.  The 


parallel  integrations  to  find  velocities  and  accelerations  may  be  avoided,  by  using 
the  following  relations,  some  of  which  are  exact,  while  others  are  only  good 
approximations . 


At  a  front  of  discontinuities,  i.e.  a  P-  or  S-front,  the  accelerations  are  of 
course  infinite,  but  the  changes  in  velocity  are  given  -  exactly  -  in  terms  of  the 
respective  stress  discontinuities, 


at  «p  =  cpp 
at  cp  =  cpg 


4iN 

U  I 

p  V  sin  q?. 

At 

p  V  sin  ep, 

(V-l) 

(V-2) 


where  the  subscripts  N,  T  indicate  normal  and  tangential  directions,  respectively. 
The  value  of  Ac^  at  a  P-front  can  be  taken  directly  from  the  numerical  computations 
for  the  stresses.  The  value  At  can  easily  be  confuted  and  is  given  in  the  tables 
listing  numerical  results. 


In  con-'  'nuous  elastic  regions  velocities  do  not  change,  while  accelerations 
vanish.  Inelastic  regions  being  reasonably  narrow,  one  may  disregard  tangential 
accelerations  and  changes  in  velocity,  while  the  normal  acceleration  may  be  assumed 
to  be  uniform  in  the  region,  giving  a  linear  change  in  velocity.  The  total  change 
in  velocity,  Au^  ,  in  an  inelastic  region  of  extent  A«p,  may  be  found  from  the  change, 
A a,  in  the  principal  stress  g ^  at  both  ends  of  the  region 


where  cp^  is  an  angle 
the  values  <p^  at  the 


Au 

a/ 

A  a 

N 

p  V  sin  sp 

B 1 

defining  the  location  of  the  plastic  region, 
end  points.  The  acceleration  is 


/V.5^ 

\  *  f 


say  the  mean  of 
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v 

Except  for  values  of  —  very  close  to  unity,  the  plastic  regions  are  narrow, 

cp 

so  that  the  above  approximations  are  quite  satisfactory  when  determining  shock 
factors . 

c)  Conclusions. 

A  method  has  been  presented  and  numerical  results  have  been  tabulated  for  the 
problem  of  a  hal_ -space  of  a  von  Mises  material  subject  to  a  step  load  progressing 
with  superseismic  velocity,  V  >  Cp  . 

With  regard  to  application  of  the  results  to  protective  construction,  it 
must  be  emphasized  that  the  effects  of  a  step  load  and  those  of  a  decaying 
pressure  differ  necessarily.  The  results  for  the  step  load  may  be  used 
approximately  for  a  decaying  load  only  for  points  in  the  vicinity  of  the  front 
of  the  surface  load  and  soon  after  passing  of  this  front.  The  solution  for  the 
step  load  approximates  the  solution  for  the  decaying  load  reasonably  within  a 
distance  D  from  the  front,  Pig.  14.  This  distance  is  the  distance  in  which  the 
peak  pressure  decays  by  about  10-20$. 

It  should  be  noted  that  the  present  report  emphasizes  numerical  results. 

More  theoretical  and  academic  matters,  e.g.  concerning  uniqueness,  are  treated 
in  [5]. 


1 
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APPENDIX  A  -  Approximate  Forms  of  Equations  (111-18-21  and  27 )» 


Using  the  relations 


0  *  3  +  A  (A-l) 

<?  -  9  +  e 


A  set  of  approximate  relations  will  be  derived  to  replace  Eqs.  (111-18-21  and  2?) 
in  the  vicinity  of  «p  -*  where  the  numerical  integration  of  the  above  equations 
encounters  difficulties,  A  and  c  are  to  be  considered  small,  but,  to  obtain 
sufficient  range  of  validity  of  this  approximation  the  first  two  significant 
terms  in  each  variable  are  retained. 


Replacing  «  -  w  -  cp  by  the  more  convenient  variable 


5  =  x  -  x 


(A-2) 


where  X,  X  are  defined  by  Eqs.  (11-34),  (HI-33)>  respectively.  To  the  order 
of  the  approximation  used  one  finds 


*  2(1+V) 

5  -  “  30^1 


x  .  3^(1-*)  , 

...X  -2  X  «  2  ,,  .  .  * 


(1+V)  c 


2Cp  (1+V) 


Similarly,  approximate  expressions  for  cos  2y,  sin  2y  are 

cos  2y  *  -  1  +  27p 
sin  2y  »  -  27} 


(A-3) 


(A-4) 


Introducing  Eqs.  (A-l)  to  (A-4)  into  the  determinantal  equation  (III -13 )  and  the 

t 

differential  equation  (111-20)  for  9  ,  retaining  only  appropriate  terms  yields 


*1^ 


+  a( 


■s 


5 


(A-5) 
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S3 

d« 


«  1  + 


% 

*5 


H  L 


(a~6  ) 


Combination  of  Eqs.  (III-17,  18  and  19)  permits  formulation  of  a  relation  for 


dA 

d« 


(A-7) 


where 


^  -  4#  -  i  (1-5V)  E  +  4(l-2v)  Is  (a-8) 

<2  -  -  32(l+'<)  -  2U(1-2V)  E  -  |  (17+8v)  A  -  8(l-2v)  Sa  (A-9) 

»3  -  8(l-8v)  -  48(l-2v)  5  +  Mi-8'0  A  -  24(l-2v)  5a  (A-10) 

-  Ml+v)  +  |  A  +  3(l-2v)  5  -  97p  -  3lfA  (A-ll) 

«  2(1-8v)  -  3(5-3U6v)  5  +  A  -  (5-l6v)  (A  + 

+  18  S2(1-2v)  +  6(1-2v)  52A  (A-12) 

«  12  +  6A  +  d2  (A-13) 

a?  *  3(1-2V)  5  +  d  +  2(l+v)  1?  +  (1-2V)  5A  (A-l4) 

ag  -  -  4(l+v)  -  |  A  +  9^2  -  3(1-2V)  5  +  3  Alp  (A-15) 


Equations  (A-5  to  7 )  govern  the  solution  in  terms  of  the  three  unknowns  A  and 
L  >  0. 


The  derivative  of  becomes 


"i  3  U+v) 

d«  1-2V 


1 


2  *9 

^  ®10 


(»1  +  ®2)  L 


(A-16) 


where 


-  3^  - 


(A-17) 
(A-18) 

The  quantities  s^  and  can  be  found  from  the  yield  condition  and  from  the 
relation  0  =  3  +  A,  once  A  is  founds 

It  may  be  shown  easily  that  Fqs.  (A-5,6)  reduce  to  their  counterparts 
presented  in  [$]  when  §,  A  and  T]  become  small  enough  to  permit  retention  of 
leading  terms  only.  Equation  (A-7)  also  reduces  to  its  counterpart  in  [5] 
providing  the  relationship  between  1)  and  A  given  by  the  simplified  determinant al 
equation  {Eq.  (3-35)  of  [5]}  is  employed. 


»  9(l-2v)  §  +  (1+V)  A  +  6(l+v)  if  +  A2  +  6(l-2v)  AS  + 

+  (l-2v)  A2S  +  2(l+v)  Alf 

aiQ  =  -  „  2(1-2V)  5  +  2Alf 
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APPENDIX  B  -  Solutions  for  v  =  g  near  cp  a  cp^  a  jp. 

1 

For  v  *  g  the  asymptotic  solutions  valid  for  cp  <  ^  *  cp  can  be  obtained  from 
a  simplification  of  the  differential  equations  (A-5,  6  and  7)  where,  after  sub¬ 
stitution  cf  v  «  g  ,  only  the  leading  terms  have  been  retained: 


36  Tf  +  £  5A2  =  365s 

(B-l) 

A'  »  -  (65  +  A)  L 

(B-2 ) 

* =1"!l 

(B-3) 

Using  tnis  asymptotic  solution  for  <p  <  cp  ■  cp  and  a  double  accuracy 

'  o 

integration  of  the  original  differential  equations  (III-18  to  27)  for  cp  >  q>  ■  qj, 

b 

solutions  could  be  constructed  successfully.  As  a  typical  case,  Figs.  15  and  l6 
show  the  values  of  A  and  I]  as  functions  of  c  for  V  s  1.25  Cp  .  There  are  only 
two  extremely  narrow  neutral  regions  between  the  two  plastic  regions  .and  the 
S-front.  Further,  the  discontinuity  in  shear  at  the  S-front  becomes  extremely 
small,  approximately  zero. 
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TABLE  1  (For  definition  of  regions,  see  Fig.  13) 


v  =  O.o  V  =  1.25  cp  cpp  =  126.87°  <ps  =  145.55°  9  -  152.49° 


pe 

p  *1.4l'k 

p_  <  p  <  p_ 

*0 

p  *2.17  k 

pL 

p  «2.40'k 

0 

*  po  <  % 

p  »  3*20  k 

0 

pH 

po  *8.02  k 

SSI 

-  0.876 

-  1.58 

-  1.73 

-  1-73 

-  1.73 

Og/k 

0 

0 

0 

0 

0 

Vk 

-  O.876 

-  1.58 

-  1.73 

- 1.73 

-  1.73 

3 

3.00 

3.00 

3.00 

3.0C 

e 

36.87° 

36.87° 

36.87° 

36.87° 

36.87C 

*1 

Hot  Applicable 

Not  Applicable 

140.29° 

140.29° 

140.29° 

’a 

If 

It 

140,29° 

143.51° 

144.08° 

ox/k 

ft 

ft 

Not  Applicable 

-  1.86 

-  1.89 

V'k 

It 

It 

It 

■H 

-  0.211 

Vk 

It 

tt 

If 

-  2.11 

-  2.22 

3 

If 

If 

tf 

huh 

e 

II 

" 

ft 

46.29° 

48.6° 

4T4 

-  1.17 

-  1.11 

-  1.05 

-  0.536 

-  o.4o4 

Vk 

m— — — — 

Hot  Applicable 

-  1.68 

-  1.73 

-  1.86 

-  1.89 

Vk 

It 

0.096 

0 

-  0.170 

-  0.211 

Vk 

It 

-  1.58 

-  1.73 

-  2.11 

-  2.22 

3 

II 

3-37 

3.00 

2,72 

2.69 

0 

It 

. 

78.00° 

74.23° 

64.81° 

62.5° 

^3 

161.39° 

156.8° 

152.42° 

O 

161.39° 

161.39° 

161.39° 

l6l. 40° 

161. 4o° 

Vk 

-  1.41 

-  2.17 

-  2.40 

-  3.20 

-  8.02 

Vk 

0.529 

-  0.333 

-  1.36 

-  6.24 

Vk 

-  0.876 

-  3*08 

-  3.77 

-  6.16 

-  20.7 

3 

6.62 

4,06 

3.98 

4.01 

4.o4 

B 

90° 

90° 

90° 

90° 

90° 

Region  D  j  Region  C  j  j  Region  B  !  Region 


TABLE  2  (For  definition  of  regions,  see  Fig.  13 ) 
v  =  0.0  V  =  1.5  tpn  =  130.19°  =  151.87°  $  =  157-36 


?jA 


°2/k 


Vk 


-  1.18 


-  1.44 


-  1.73 


pT  <  p  < 

■■I. 


-  1.73 


-  1.18 


-  1.44 


3.00 


-  1-73 


•3.00 


48.19°  48.19°  48.19°  - 


ot  Applicableffot  Applicable  1^.90° 


148.90° 


o2/k 


-  1.05 


-  0.978 


-  0.796 


*Vk  tot  Applicable!  .  -  I.63  -  1.73 


94 


Vk 


163.48° 


163.48° 


-  1.53 


0.352 


-  1.18 


161.72° 


163.530 


158. 88c 


-  0.135 


-  2.11 


-  2.38 


-  0.571 


-  3.70 


-  1*73 


3.00 


48.19° 


148.90° 


150.02° 


ot  Applicable  “  1*78 


-  0.0689 


-  1.89 


2.88 


52.15° 


-  O.569 


-'I.78 


157.89 


163.48° 


-  2.71 


-  0.893 


-  4.67 


-  1.73 


-  1.73 


3.00 


48.19° 


148.90° 


151.09° 


“  1.83. 


-  0.124 


-  2.02 


2.81 


56.38° 


-  0.327 


-  1.83 


-  0.124 


-  2.02 


2.81 


67.37° 


157.34° 


163.48° 


-  6.57 


-  4.79 


-  16.32 


-  39 


TABLE  3  (For  definition  of  regions,  see  Fig.  13) 


v  =  0.0  V  =  2.00  cp  cpp  =  150.00 


<ps  =  159.30 


1.65s 


<Vk 

Vk 


- 1.44 


- 1.44 
3.00 
6o.ooc 


•Q„  <  t)  <  p„ 

-c 

P0  «  1.954 
-  1.62 


-  1.62 
3.00 
60.00° 


■L 

2.32k 


-  1.73 


-  1.73 
3.00 
60.00° 


tot  Applicable  Not  Applicable  157*95" 
"  "  157. 95C 


pL  <  po  <  p 
PQ  «  2.7& 

-  1.73 


-  1.  /. 


3.00 

60.00° 

157*95° 

158.66° 


I 


ajk 

JL 

c2/k 

Jx/k 


-  0.831 


0.696 


Og/k 


tot  Applicable  -  I.69 


tot  Applicable  -  1.77 

"  -  0.0495 

"  -  I.85 

" _ 2.92 

" _ 63.57° 

-  0.552  -  0.342 

-  1.73  -  '1-77 


-  0.0495 


1.73 


3-00 


60.00° 


157-95° 


158.95° 


-  1.78 


-  0.0671 


-  1.89 


2.90 


65.40° 


-  0.237 


-  1.78 


0.0671 


73-26° 


163.22° 


166.54° 


-  6.61 
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TABIiE  4  (For  definition  of  regions,  see  Fig.  13) 


V  =  0.0  V  =  4.00  cu  cp_  =  165.52°  oc  =  169.82°  cp  «.  171.70° 

ir  ir  o 


PE 

P„  -  1.73* 

PE  <  po  <  PL 
Pft  *  1.91k 

p0  *  2.17k 

pL  <  p0  4  PH 

po  -  a. 60k 

pn  -  7.12k 

II 

turn 

'  -  1.66 

-  1.71 

-  1.73 

- 1.73 

-  1.73 

< 

Vk 

0 

0 

0 

0 

_ SL _ _ 

fl 

0 

■a 

Vk 

-  1 .66 

- 1.71 

-  1.73 

< 

- 1.73 

-  1.73 

$ 

3 

3.00 

3.00 

3*00 

3.00 

3.00 

1 

e 

75*52° 

75.52° 

75.52°- 

75.52° 

75.52° 

11 

*l 

[Jot  Applicable 

Not  Applicable 

169.61° 

169.61° 

■ 

*2 

tf 

ft 

169.61° 

169.74° 

169.77° 

Vk 

ft 

ft 

Not  Applicable 

-  1.74 

-  1.75. 

cq 

CTg/k 

ft 

ft 

ft 

-  0.0145 

-  0.018 

CS 

0 

■a 

Vk 

ft 

ft 

ft 

- 1.77 

-  1.77 

£ 

3 

ft 

" 

It 

2.98 

2.97 

e 

ft 

tf 

ft 

77.34° 

77*89° 

■ 

At4 

-  0.429 

-  0.341 

-  0.259 

-  0.149 

-  0.116 

ai/k 

lot  Applicable 

-  1.72 

-  1-73 

-  '1.74 

-  1.75 

0 

Vk 

ft 

0.0171 

0 

MKgMMm 

119 

-  0.18 

c 

0 

•a 

Vk 

t? 

-  I.71 

-  1-73 

-  1.77 

-  1.77 

<i) 

3 

ft 

3.06 

3.00 

2.98 

2.97 

L 

e 

‘ft 

66.93° 

84.11° 

82.29° 

81.74° 

T3 

172.38° 

172.02° 

171.78° 

171.71° 

171.70° 

* 

*4 

172.38° 

172.38° 

172.38° 

172.38° 

172.38° 

oL/k 

-  I.71 

-  1.91 

-  2.17 

-  2.60 

-  7.12 

« 

a 

0 

*9) 

<?g/k 

0.0455 

-  0.157 

-  0.425 

-  0.856 

-  5*37 

Vk 

-  1,66 

-  2.26 

— 

-  4.35 

-  17.9 

<D 

PS 

3 

3.16 

3.13 

•  3.11 

3;il 

3.12 

9 

90° 

90° 

90° 

. . 

90° 

'  90° 

TABLE  5  (For  definition  of  regions,  see  Fig.  13) 


1 

1 

'! 


V  =  0.125  v  =  1*25  Cp  9p  =  126.87°  cps  =  148.42°  cp  =  148.42° 


P„  1 
Pn  *  1.721c 

-JE  <  ?o  <  Sl 
P„  =  2.09k 

PL 

p  =  2.36lc 

0 

»L  <  »o  <  % 
PQ  *  3.651C 

% 

p  *  6.l4lc 

! 

ox/k 

-  1.36 

- 1.75 

-  2.02 

-  2.02 

-  2.02 

< 

a2/k 

-  0.194 

-  0.2',! 

-  0.289 

-  0.289 

-  0.289 

G 

0 

fcD 

Vk 

-  1.75 

-  2.26 

-  2.60 

-  2.60 

-  2.60 

<D 

K 

3 

1 

3.00  i 

3.00 

•3.00 

3-00 

3.00 

0 

36.87° 

36.87° 

36.87° 

36.87° 

36.87° 

*1  : 

tot  Applicable 

tot  Applicable 

140.38° 

140.38° 

140.38° 

V2 

It 

It 

140.38° 

148.32° 

148.42° 

Vk 

ft 

II 

tot  Apj.  ,.i cable 

-  2.64 

-  5.15 

PQ 

CT2/k 

If 

» 

IT 

-  0.958 

-  3.47 

0 

u 

v* 

II 

II 

II 

-  4.47 

-  12.08 

d; 

« 

3 

II 

It 

If 

2.72 

2.99 

1 

i 

i 

0 

If 

If 

II 

58.19° 

58.42° 

4Ti 

-  1.24 

-  1.23 

-  1.16 

-  O.0135 

0 

Vk 

tot  Applicable 

-  I.90 

-  2.02 

-  '2.64 

-  5.15 

o 

Vk 

II 

-  0.102 

-  0.289 

-  0.958 

-  3M 

c 

o 

u 

Vk 

If 

-  2.26 

-  2.60 

-  4.47 

~  12.08 

<u 

0 

IV 

3.59 

3.00 

2.72 

2.99 

• 

e 

II 

84.65° 

79-97° 

58.64° 

58.1*2° 

■ 

*3 

160.90° 

158.62° 

156.59° 

148.49° 

148,42° 

l60.900 

160.83° 

160.80° 

160.85° 

160.85° 

°iA 

-  1.72 

-  2.09 

-  2.36 

-  3.65 

-  6.14 

« 

Vk 

O.I65 

-  0.248 

-  0.541 

-  1.80 

-  4.34 

n 

o 

n> 

Vk 

-  1*75 

-  2.83 

-  3.65 

-  7.51 

-  15^1 

« 

3 

4.85 

4.06 

3.85 

4.14 

4.25 

e 

90° 

90° 

90° 

90° 

90° 

Region  D  '  Region  C  Region  B  Region 
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TABLE  6  (For  definition  of  regions,  see  Fig.  13) 
v  =  0.125  V  =  1.50  cp  cpp  =  138.19°  cps  =  154.12°  cp  =  154.12° 


V 

PE 

P  a  1.82k 

O 

»E  <  *0  <  *L 

p  =  2.04k 

0 

PL 

Pn  *  2.35k 

*L  <  Po  < 
p0  -  3.26k 

pH 

Pn  *  5* 95k 

°iA 

-  1.57 

-  1.77 

-  2.02 

-  2.02 

-  2.02 

Vk 

-  0.224 

-  0.253  ' 

-  0.289 

-  0.289 

-  0.289 

Vk 

-  2.02 

-  2.28 

-  2.60 

-  2.60 

-  2.60 

3 

3-00 

3.00 

■3‘00 

3.00 

3.00 

e 

48.19° 

48.19° 

48.19°  • 

48.19° 

48.19° 

■Jot  Applicable 

'Jot  Applicable 

149.05° 

149.05° 

149.05° 

*2 

It 

tt 

149.05° 

153.81° 

154.12° 

) 

Vk 

ft 

"  Lot  Applicable 

-  2.46 

-  5*14 

CT2/k 

If 

ft 

tt 

-  0.759 

-  3*44 

Vk 

If 

ft 

tt 

-  3.92 

-  12.01 

3 

1! 

Jf 

it 

2.80 

2.99 

e 

It 

tt 

63=,  16° 

64.12° 

K 

-  1.08 

-  1-03 

-  0.914 

-  0.0569 

0 

> 

Vk 

. 

tot  Applicable 

-  1.91 

-  2.02 

-  "2.46 

-  5.14 

Vk 

ft 

-  0.113 

-  0.289 

-  0.759 

-  3^44 

JXA 

tf 

-  2.28 

-  2.60 

-  3.92 

-  12.01 

e 

It 

3.55 

3.00 

2.80 

2.99 

e 

tf 

86.22° 

80.06° 

65.08° 

64.12° 

CP- 

0 

•  /■  — 

LO'd.O^ 

-r>  O 

xox . yo 

-/-O 

±5*4.  ^0 

. -O 

*4 

162.89° 

162.940 

l62.89C 

162.93° 

162.95° 

b 

1 

Vk 

-  1.82 

-  2.04 

-  2.35 

-  3.26 

-  5-95 

OgA 

0.0259 

-  0.218 

-  0.546 

-  1.44 

-  4.16 

Vk 

-  2.02 

-  2,67 

-  3*61 

-  6.32 

-  14.47 

3 

4.12 

3.83 

3*66 

3<8'o 

3.89 

e 

90° 

90° 

•90° 

90° 

0 

90 

Region  D  '  Region  C  Region  B  Region 
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TABLE  7  (For  definition  of  regions,  see  Fig.  13 ) 
v  =  0.125  V  -  2.00  Cp  cpp  =  150.00°  cp3  =  I6O.890  cp  =  160.89° 


vk 


a2/k 


Jl/k 


0.25^ 


-  0.289 


-  2.48 


60.00 


158.16° 


-  2.6c 


60.00° 


158.16° 


159.80° 


-  2.21 


-  2.60 


60.00 


160.89° 


-  5.16 


Region  D  '  Region  C  Region  B  Region 
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TABLE  8  (For  definition  of  regions,  see  Fig.  13 ) 


v  =  0.125  v  =  4.00  cp  cpp  =  165.52 


PB  PE<Po<PL 

p  =  2.00k  p  «  2.09k 

*  A  I  •  r*\ 


2.16k 


vk 

-  I.96 

~  2.00 

-  2.02 

c2/k 

-  0.280 

-  0.286  * 

-  0.289 

Vk 

-  2.52 

!  -  2.57  * 

-  2.60 

0 

3-00 

3.00  j 

!  3.00 

0 

75° 52° 

j 

75.52° 

I  75.52° 

*1 

tot  Applicable 

tot  Applicable 

t 

169.87° 

*2 

1 

1  U 

ft 

169.87° 

Vk 

Vk 

Vk 

e 

~ e 

At4 

giA 

Vk 

Vk 

s 

0 


tot  Applicable 


-  0.430 
tot  Applicable 


-  O.36O 

-  2.01 

-  0.273 

-  2.57 
3.04 

87.55° 


-  0.304 

-  2.02 
-  0.289 
-  2.60 

3.00 

85.64° 


*3 

171.56° 

171.59° 

171.31° 

*4 

171.96° 

171.96° 

171.96° 

Vk 

-  2.00 

-  2.09 

-  2.16 

tfgA 

-  0.246 

-  0.341 

-  o.4n 

Vk 

-  2.52 

-  2.80 

-  3.01 

0 

3.12 

3.H 

3.11 

9 

90° 

90° 

. 

VO 

O 

O 

cps  =  170.58 

PL  <  po  <  p 
pQ  *  2.30k 

-  2.02 

-  0.289 

r 

-  2.60 
3.00 

75*52° 

169.87° 

170.20° 

-  2.09 

-  0.364 

-  2.81 

-  2.97 
77*92° 

-  0.160 
-'2.09 

-  0.364 

-  2.81 

2*97 

83.24° 

170.96° 

171.96° 

-  2.30 

-  0*555 

-  3*44 
3*11 


-  2.02 


-  0 


-  2.60 


169.87° 


170.58° 


-  5.16 


-  3*42 


-  12.00 


3*00 


80.60° 


0 


-  5.16 


-  3.42 


-  12.00 


3*00 


80.60° 


0.58° 


171.960 


-  5.44 


-  3.70 


-  12.87 


3.14 
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TABLE  9  (For  definition  of  regions,  see  Fig.  13) 


v  =  0.25  V  =  1.25  =  126.87°  <pc  =  152.49°  <p  »  143. 40° 

*  4  O 


1 

PE 

P0  -2.50  k 

'  ! 

p„  ^  0  <  p_ 

Pc  -2.61  k 

PL 

Pn  “2.65  k 

H  <  p0  <  % 

PD  -3-84  k 

PH 

pft  -10.5  k 

1 

Region  A 

EM 

-  2.43 

-  2.55 

-  2.60 

-  2.60 

-  2.60 

cr^/k 

d 

-  0.810 

-  0.850 

-  0.866  -  0.866 

-  0.866 

Vk 

-  4.05 

-  4.25 

-  4.33 

-  4.33 

-  4.33 

0 

3*00 

3.00 

3.00 

3.00 

3.00 

e 

36.87° 

36.87° 

36.87° 

36.87° 

36.87° 

*1 

tot  Applicable 

Not  Applicable 

139. 39° 

139.39° 

139.39° 

*2 

If 

tt 

139. 39° 

142.51° 

143.3960 

Region  B 

Oj/k 

tl 

"  |Not  Applicable 

-  3.51 

-  10.1 

a2/k 

»r 

tt 

It 

-  1.81 

-  8.35 

Vk 

tr 

It 

•t 

-  7c08 

-  26.8 

0 

tt 

tt 

it 

2.82 

00 

• 

0 

0 

0 

it 

It 

tt 

47.^6° 

53.4° 

H 

At4 

~  1.36 

-  1.35 

-  1.35 

-  0.851 

-  0.550 

Region  C 

Vk 

tot  Applicable 

-  2.57 

-  2.60 

-  3.51 

-  10.1 

Vk 

tt 

-  O.83O 

-  0.866 

-  1.81 

-  8.35 

Vk 

tt 

-  4.25 

-  4.33 

-  7.08 

-  26.8 

0 

tl 

3.07 

3.00 

2.82 

3.00 

8 

tt 

88.60° 

88.11° 

77.52° 

71.6° 

co_ 

b 

.  f  - J  0 

iOv  (Ox 

1  /f  r\  /■\r*Q 

JLDV.  V? 

156.12° 

1  ell  l.rt° 

JLJH*  Hv 

’4 

160.61° 

160.61° 

160.59° 

160.88° 

161.01° 

Region  D 

Vk 

-  2.50 

-  2.61 

-  2.65 

-  3.84 

-  10.5 

Vk 

-  0.735 

-  0.851 

-  0.897 

-  2.03 

-  8.69 

Vk 

-  4.05 

-  4.38 

-  4.50 

~  8.07 

-  28.1 

0 

3.28 

3.21 

3.19 

3.74 

4.12 

0 

QO° 

7V 

90° 

90° 

90° 

90° 

J 
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TABLE  10  (For  definition  of  regions,  see  Fig.  13) 


V  a  0.25  v  =  1.50  Cp  <pp  =  138,19°  cps  =  157.30°  cp  »  150.20° 


»E 

P0  -  2.Wk 

>E  <  *0  <  »L 
Pn  «  2.63k 

pL 

pn  »  2.69k 

V<rH 

-  4.i3Uc 

PH 

-  11.6k 

Region  A 

°l/k 

“  2.39 

-  2.54 

-  2.60 

-  2.60 

-  2.60 

Vk 

-  0.797 

-  0.846 

-  0.866 

-  0.866 

-  0.866 

Vk 

-  3.98 

-  4.23 

-  4.33 

-  4.33 

-  4.33 

0 

3.00 

3.00 

3.00 

3.00 

3-00 

e 

48.19° 

48.19° 

48.19° 

48.19° 

48.19° 

Mot  Applicable 

Not  Applicable 

148. 06° 

148. 06° 

148.06° 

*2 

ft 

It 

148.06° 

150.13° 

150.20° 

Region  B 

Vk 

If 

It 

Not  Applicable 

-  3-99 

-  11.25 

a2/k 

II 

If 

II 

-  2.27 

-  9.53 

JXA 

II 

ft 

II 

-  8.51 

-  30.31 

3 

II 

If 

ft 

_ 2.91 

3.00 

e 

1! 

II 

If 

58.51° 

60.21° 

□ 

At4 

-  1.13 

-  1.09 

-  1.07 

"  0.523 

-  0.425 

Region  C 

Vk 

Mot  Applicable 

-  2.56 

-  2.60 

‘  3.99 

-  11.25 

Vk 

If 

-  0.820 

-  0.866 

-  2.27 

-  9.53 

J3/k 

II 

-  4.23 

-  4.33 

-  8.51 

-  30.31 

0 

It 

3.09 

3.00 

2.91 

3c00 

6 

If 

87.61° 

86.54° 

76.21° 

74.52° 

*3 

162.95° 

162.14° 

■ 

161.78° 

158.83° 

158. 1(3° 

fk 

162.950 

162.95° 

162.95° 

163.09° 

163.130 

Region  D 

Vk 

-  2.48 

~  2.63 

-  2.69 

-•  4.31 

-  11.60 

V* 

~  O.705 

-  0.861 

-  0.925 

-  2.50 

~  9.78 

Jl/k 

-  3- 98 

-  4.43 

-  4.6l 

-  9.47 

-  30.00 

0 

3.34 

3-28 

3.26 

3.66 

3.78 

0 

90° 

909 

90° 

90° 

90° 

Region  D  I  j  Region  C  Region  B  I  Region 


TABLE  11  (For  definition  cf  regions,  see  Fig.  13) 

V  =  0.25  V  =  2.00  cp  cpp  =  150.00°  cps  =  163.220  <p  a  158.12° 


PE 

Pn  -  2.52k 

tE<I’o<  pL 
P..  ®  2.6lk 

PL 

P^  »  2.68k 

»L  <  ®o  <  *H 
?e  "  4  Aik 

PH 

p.  ■  11.4k 

BKa 

-  2.1+6 

-  2.54 

-  2.60 

-  2.60 

-  2.60 

-  0.820 

-  0.84 

-  0.866 

-  0.866 

-  0.866 

V* 

-  4.10 

-  4.23 

-  4.33 

-  4.33 

-  4.33 

a 

3.00 

3.00 

3.00 

3.00 

3.00 

e 

60.00° 

_ 

60.00° 

60.00° 

60,00° 

60.00° 

Not  Applicable 

157.22° 

157.22° 

157.22° 

*2 

ft 

If 

157*22° 

i58.ll0 

158.12° 

If 

It 

Not  Applicable 

-  4.20 

-  11.16 

CT2/k 

!» 

If 

It 

mm 

-  9.43 

JiA 

»» 

If 

It 

-  9.12 

-  30,02 

a 

If 

If 

It 

-  2.96 

3.00 

e 

tl 

It 

II 

67.46° 

68.12° 

h 

-  0.852 

-  0.811 

-  0.771 

-  0.672 

-  0.306 

Vk 

Not  Applicable 

-  2.56 

-  2.60 

-  4.20 

-  11.16 

Og/k 

If 

-  0.818 

-  0.866 

-  2.47 

-  9.43 

Vk 

If 

HHfli 

-  4.33  ■ 

-  9.12 

-  30.02 

a 

If 

3.10 

3-00 

2.96 

3.00 

e 

If 

88.17° 

86.44° 

| 

78.32° 

<Po 

mLsaxtmm 

-*-w  S* 

*70° 

16-3  £o° 

166.21° 

166.21° 

166.21° 

166.26° 

166.27° 

V‘ 

-  2.52  • 

•>  2.61 

-  2.68 

-  4.4l 

-  11.38 

as/k 

-  0.759 

-  0.848 

-  0.923 

-  2.63 

-  9.60 

Vk 

-  4.10 

-  4.36 

-  4.58 

-  9.77 

-  30.69 

a 

mam 

■m 

■m 

3.40 

3.44 

0 

_ _ _ 

90° 

90^ 

90° 

90° 

90° 

Region  B  Region  C  Region  B  Region 
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TABLE  12  (For  definition  of  regions,  see  Fig.  13) 

V  a  0.25  v  a  4.00  cp  <p  =  165.52°  «PS  =  171*70°  $  =■  169.26° 


PE 

Pn  «  2.58k 

PE  <  po  <  PL 
P0  «  2.62k 

pL 

p  »  2.63k 

0 

PL<P0<PH 
PQ  *  4.36k 

pH 

Pft  -  11.3k 

Region  A 

Km 

1 

"  2.59 

-  2.60 

-  2.60 

-  2.60 

Cfg/k 

-  0.854 

-  0.863 

-  0.866 

-  0.866 

-  0.866 

JiA 

-  4.27 

-  4*33 

-  4.33 

3 

3*00 

3.00 

3*00 

3*00 

3.00 

e 

75*52° 

75*52° 

75*52° 

75*52° 

75.52° 

1 

*1 

tot  Applicable 

tot  Applicable 

169*14° 

169.14° 

169.14° 

*2 

tl 

If 

.  _ 

l69*l4° 

169.26° 

169.26° 

Region  B 

°l/k 

it 

If 

tot  Applicable 

-  4.29 

-  11.22 

°SA 

It 

If 

If 

-  2.56 

-  9*49 

Vk 

t» 

If 

If 

-  9*4o 

-  30.2 

3 

tt 

If 

«f 

mam 

waam 

e 

If 

If 

If 

msm 

■ 

4T4 

-  0.431 

-  O.388 

-  0.371 

-  0.337 

-  0.147 

Region  C 

Vk 

Not  Applicable 

>  -  2.59 

-  2.60 

-  4.29 

-  n.22 

Cg/k 

If 

-  0.859 

-  0.866 

-  2.56 

-  9.49 

Vk 

M 

-  4.31 

-  4.33 

-  9.40 

-  30.2  . 

3 

If 

3.01 

3*00 

2.99 

3.00 

e 

If 

■■ 

11111121 

84.35° 

84.14° 

Fa 

1 

172.11° 

172.05° 

171.78° 

171.76° 

El 

172. 27° 

172.27° 

172.27° 

172.27° 

172.27° 

p 

0 

% 

& 

cr-j/k 

-  2.58 

-  2.62 

-  2.63 

-  4.36 

-  11.3 

Vk 

■ni 

-  0.877 

-  0.891 

-  2.61 

-  9.55 

Vk 

msm 

-  4.43 

-  9.61 

-  30.42 

9 

3.06 

3.06 

3*06 

3.10 

3.11 

9 

90° 

90° 

90° 

90° 

90° 
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TABLE  13  (For  definition  of  regions,  see  Fig.  13) 


v  =  0.35  v  =  1.25  Cp  cpp  =  126.87°  cps  =  157*40°  $  =,  138.27° 


rE 

P0  »  3* 48k 

PE  <  P0  *  pL 
P„  •  4.43k 

pL 

p  m  6 . 08k 

-  o_ 

PL  <  Po  <  PH 
PQ  -  6.72k 

PH 

p  «  11.4k 

0 

Region  A 

vk 

-  3.75 

-  3*75 

-  3.75 

~  3.75 

-  3-75 

v'k 

-  2.02 

-  2.02 

-  2.02 

BBBHI 

-  2.02 

Vk 

-  7*79 

-  7.79 

-  7.79 

-  7.79 

-  7.79 

3 

mam 

3-00 

3.00 

3.00 

3.00 

9 

36.37° 

36.87° 

36.87° 

36.87° 

36.87° 

1 

*1 

BPSSM 

136.870 

m 

136.87° 

136.87° 

*2 

136.87° 

138.10° 

138.18° 

138.27° 

M 

C 

1  O 

!  *4 

Vk 

Not  Applicable 

-  4.57 

~  6.08 

-  6.69 

-  11.31 

°SA 

II 

-  2.84 

-  4.35 

-  4.97 

-  9.58 

JjA 

ft 

-  10.23 

-  14.76 

-  16.61 

»  30.48 

1  <y 

1  « 

» 

1 

3 

ft 

2.94 

2.93 

2.94 

2.99 

e 

II 

42.49° 

44.80° 

45.79° 

48.09° 

1 

1 _ 

At4 

-  1.18 

-  1.18 

- 1.22 

- 1.18 

-  1.08 

Region  C 

_  _ 

Vk 

tot  Applicable 

tot  Applicable 

-  6.69 

-  11.31 

yk 

H 

II 

it 

-  4.97 

-  9.58 

Vk 

SI 

II 

!• 

-  16.61 

-  30.48 

3 

If 

II 

H 

2.94 

2.99 

e 

II 

ft 

it 

89.01° 

86.71° 

■ 

T3 

11 

If 

-i/r-i  c-,® 
JLOJu.U  1 

i6"l 

*4 

te 

II 

161.97° 

162.00° 

162.11° 

O 

fl 

0 

axA 

-  3*48 

-  4.43 

-  6.08 

-  6,72 

-  11.39 

-  2.29 

-  2.98 

-  4.35 

-  4.98 

-  9-63 

Vk 

-  7.79 

-  10.23 

-  14.76 

-  16.68 

-  30.71 

3 

2.06 

msm 

2.93 

3.02 

3.24 

6 

90° 

90° 

90° 

j 

90° 

j  90° 

Region  D  j)  Region  C  Region  B  Region 


i 
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TABLE  14  (For  definition  of  regions,  see  Fig.  13) 


V  a  0.35  V  =  1.50  Cp  <Pp  *  138*19°  <pg  =  181*32°  <P  a  II46.3I 


i 


t 

t 


agy»a«*y 


Region  D  j  j  Region  C  }  }  Region  B  j  Region 


t 
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TABLE  15  (For  definition  of  regions,  see  Fig.  13) 

v  =  0.35  v  =  2.00  cp  cpp  =  150.00°  cps  =  166.10°  cp  »  155.42° 


PE 

Pn  -  3.72k 

?E  <  po  <  PL 
pn  «•  4.34k 

PL 

PQ  -  4.94k 

PL  <  Po  *  PH 
PQ  «  5.7HC 

PH 

■  11.8k 

_ LO 

mm 

-  3-75 

-  3.75 

-  3.75 

"  3.75 

WBBBSKL 

Vk 

-  2.02 

-  2.02 

-  2.02 

-  2.02 

1 

Vk 

-  7.79 

-  7.79 

-  7.79 

-  7.79 

-  7.79 

3 

3.00 

3.00 

3-00 

3.00 

3.00  ■ 

8 

6o° 

6o° 

60° 

6o° 

6o° 

*1 

155.21° 

155.21° 

155.21° 

155.21° 

155.21° 

*2 

155.21° 

155.30° 

.  -  -  -  .  - 

155. 35° 

155-39° 

155.42° 

-Vk 

Not  Applicable 

-  4.36 

-  4.94 

-  5.69 

-  11.75 

<*2/k 

It 

-  2.63 

-  3-20 

-  3.96 

-  10.02 

Vk 

ft 

-  9.60 

-  11.37 

-  13.60 

-  31.78 

3 

If 

2.99 

.  .  .*•»  . 

2.98 

3.00 

3 

II 

61.52° 

62.59° 

63.58° 

65-35° 

m 

-  0.848 

-  0.820 

-  0.801 

-  0.732 

-  0.634 

Vk 

tot  Applicable 

tot  Applicable 

Not  Applicable 

-  5.69 

-  11.75 

Vk 

It 

H 

ft 

-  3.96 

-  10.02 

Vk 

VI 

If 

II 

-  13.60 

B 

3 

II 

II 

II 

mm 

3.00 

e 

II 

II 

II 

whsm 

HH 

*3 

*1 

M 

167.56° 

167.33° 

167.05° 

% 

II 

167.58° 

167.58° 

167.59° 

°l/k 

-  3.72 

-  4.34 

-  4.94 

~  5.71 

-  n.80 

Vk 
_ _ 

-  2.06 

-  2.64 

-  3-20 

-  3.97 

-  10.05 

iH 

-  7.79 

-  9.60 

-  11.37 

-  13.67 

-  31 ->93 

3 

2.87 

2.95 

2.99 

3.05 

3.14 

8 

90° 

Q0° 

90° 

.  „  .  T1  1  m  1  ill 

90° 

90° 

Region  D  Region  C  Region  B  Region 


TABLE  16  (For  definition  of  regions,  see  Fig.  13) 
v  =  0.35  V  =  4.00  Cp  cpp  •■=  165.52°  cps  =>  173.10°  9  *  167.99° 


rE 

pft  -  3.75'k 

»E  <  50  <  PL 
Pn  "  4.15k 

pL 

PG  *  5.08k 

PL<P0<»H 
p0  -  6.82k 

% 

Pn  -  llo 3k 

Region  A 

EH 

-  3.75 

"  3.75 

-  3.75 

"  3.75 

-  3.75 

ffg/k 

-  2.02 

-  2.02 

-  2.02 

-  2.02 

JiA 

-  7.79 

-  7.79 

-  7.79 

-  7.79 

3 

3.00 

3.00 

3.00 

3.00 

9 

75.52° 

75.52° 

75.52° 

75.52° 

1 

167.97° 

167.97° 

167.97° 

167.97° 

*2 

167.97° 

167.98° 

167.98° 

167.993° 

167.994° 

.<Vk 

Not  Applicable 

-  4.15 

-  5.08 

-  6.81 

-  11.25 

II 

~  2.42 

-  3.35 

-  5.08 

V* 

tl 

-  9.00 

-11.8 

-  16.97 

9 

tf 

3.00 

3.00 

3.00 

2.99 

9 

It 

WB&M 

76.84° 

77-57° 

■ 

% 

-  0.432 

-  0.418 

-  0.394 

-  O.332 

|  Region  C 

°iA 

Not  Applicable 

Not  Applicable 

Not  Applicable 

-  6.81 

It 

ft 

ft 

-  5.08 

JiA 

tf 

ft 

ft 

-  16.97 

9 

tt 

tf 

tt 

3.00 

2.99 

e 

Iff 

Iff 

»t 

88.64° 

/NA  .<  Q 

00.24 

1 

*3 

tt 

tt 

173.33° 

173.25° 

■p4 

It 

It 

173.33° 

173-33° 

Q 

§ 

% 

« 

Vk 

-  3.75 

-  4.15 

-  5.08 

-  6.82 

CgA 

-  2.03 

-  2.42 

-  3.35 

-  5.09 

-  9.52 

JiA 

-  7.79 

-  9.00 

-n.8 

-  17.01 

-  30.32 

3 

2.98 

2.99 

3.00 

3.03 

3.03 

e 

90° 

90° 

90° 

90° 

90° 

w; 


1 

1 


SBh 


»  w 


w  » 


i 

) 
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FIG  13 

Referred  to  in  Tables  1  to  1 6 

Rote  that  not  all  regions  shown  actually  occur  for  all  the  cases  of 
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